Scalar Perturbation and Stability of Ricci Dark Energy 
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The Ricci dark energy (RDE) proposed to explain the accelerating expansion of the universe 
requires its parameter a < 1, whose value will determine the behavior of RDE. In this Letter, we 
I study the scalar perturbation of RDE with and without matter in the universe, and we find that in 

, both cases, the perturbation is stable if a > 1/3, which gives a lower bound for a theoretically. 

o ■ 

^ , Observations like Type la supernovae, CMB and SDSS et al. have strongly confirmed that our universe is accelerated 
O > expanding recently. However, since ordinary matter such as stars always attract each other due to the Newton's gravity, 
the universe can only be decelerated expanding. Thus, there must be an unknown energy component living in the 
universe, and people often call it dark energy. Experiments have indicated there are mainly about 73% dark energy 
' and 27% matter components in the recent universe, but so far people still do not understand what is dark energy from 
fundamental theory. The best candidate seems the cosmology constant including the vacuum energy, but it suffers 
the fine-tuning and coincidence problems. In order to alleviate these problems, a lot of dynamic dark energy models 
^ have been built, such as quintessence, phantom, quintom models which are basically scalar field models. Another way 

■ to explain the accelerating is to modify the Einstein gravity theory like the famous /(i?) theory and DGP model etc.. 
CIh' Actually, the cosmology constant (or dark energy) problem is in essence an issue of quantum gravity , since the 
^ , density of dark energy is inevitably related to the large vacuum energy density of the quantum field theory without 

including gravity. Considering the gravity effects, there may be some regions in which the field theory can be valid. 
The holographic principle regards the black hole as the object with maximum entropy in a given region, and from 
^-H ■ statistical physics, the entropy is a extensive quantity which proportional to the volume of such region, while the black 
K*" \ hole' entropy is proportional to the area of its surface, so in the field theory, there should exists a infrared (IR) cutoff, 

■ beyond which the field theory will be failed. However, such constraint seems a little bit loose, because it includes the 
' black hole state in the field theory. To avoid the existence of such states, Cohen et al.0 suggested that in a given 
, region with length scale L, the field's energy should be bounded by the black hole's, i.e. pL^ < LM^i, where p is 



< 



O 



the total energy density within the region and Mpi = G is the Planck mass. Applying the holographic principle 
to cosmology, Li Q has proposed the holographic dark energy model, in which the energy density of dark energy is 



^\ p = 3c^MpiL~^, namely it saturates the bound. He finds that when L = Rh, which is the future event horizon, this 

. model will be consistent with observations and meanwhile solves the coincidence problem. 
^ ' Although the holographic model based on the future event horizon is successful in fitting the current data, some 
• ^ , authors asked why the current acceleration of the universe is determined by its future. Actually, the future event 
horizon is not the only choice for the holographic dark energy model. Also motivated by the holographic principle, 
H Gao, et al.Q have proposed the Ricci dark energy (RDE) model recently, in which the future event horizon area is 
replaced by the inverse of Ricci scalar, and this model is also phenomenologically viable. 

Assuming the black hole is formed by gravitation collapsing of the perturbation in the universe, the maximal black 
hole can be formed is determined by the casual connection scale Rcc given by the "Jeans" scale of the perturbations. 
For tensor perturbations, i.e. gravitational perturbations, = Max{H + 2H^,—H) for a flat universe, where 

H = a I a is the Hubble parameter, and according to the ref.[5|, only in the case of = -ff + 2i7^ it could be 
consistent with the current cosmological observations when the vacuum density appears as an independently conserved 
energy component. As we know, in flat FRW universe, the Ricci scalar is i? = 6{H + 2H^), which means the Rcc oc R 
and if one choices the casual connection scale Rcc as the IR cutoff, the Ricci dark energy model is also obtained. For 
recent progress on Ricci dark energy and holographic dark energy, see ref. [1] Hi ■ The energy density of RDE in flat 
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universe reads 



PR 



(1) 



where we have set SnG = 1 and a is a dimensionless parameter which wih determine the evolution behavior of RDE. 

In the following, we will investigate the scalar perturbation of RDE model. At first, we only consider the case when 
the universe is dominated by RDE and analytically solve the equation of motion for the perturbation. And then we 
consider the case when matter comes in, namely, the universe dominated by both RDE and matter, the equation 
of motion for perturbation is too complicated to be solved analytically, but it can be solved in the long-wave and 
short-wave limit. And we also solve it numerically. In both cases, it shows that, the perturbation is stable if a > 1/3 
and this result gives a lower bound for a theoritically. 

The linear scalar perturbation of the flat FRW metric in longitudinal (Newtonian) gauge is given by : 



ds^ = - (1 + 2$) dt"^ +a^{l- 2'if )S,jdx'dx^ 



(2) 



and we assume that the perturbation has spherical symmetry, namely, there is no anisotropic stress to the linear-order 
of the perturbation, then from the off-diagonal ij perturbed Einstein equations one can see <& = ^E". Therefore the 
perturbation of RDE is 



rv2 



$ - 3$ - 15i/$ - 6(ij + 2i?2)$ 



and the 00 component of perturbed Einstein equations is given by 



-3i7^$ - 3i/$+ ^$ 



= 5 PR + 6p„ 



(3) 



(4) 



where 5pm denotes the perturbation of matter. 

RDE without matter: First, let us consider the case when the universe is dominated by RDE only, then the 
Friedmann equation reads: 



tt2 PR 

H = — = a 



H + 2H^ = a 



2W 



where prime denotes the derivative with respect to a; = In a and hereafter we set ag = 1. 
above equation is H'^ = 7?oe~^(^~°)^ and the energy density of RDE is pn = SH^e^^'y'^^"^'' 
conservation law we get it's equation of state as 



w 



1 /2 

-3U-^ 



(5) 

The solution of the 
By using the energy 

(6) 



and to make the universe accelerated expanding (w < —1/3), it requires a < 1. Since in the linear theory of cosmology 
perturbations all Fourier modes evolve independently, we will focus on each mode labeled by its comoving wave number 
k. By using the background equation ([5|) and its solution, we obtain the following perturbation equation in momentum 
space for RDE without matter from cq.([l]) : 



3 - 



(2 



3a 



-2(^)-$, 



0. 



(7) 



Here, the second term proportional to <I>'^ is a friction term, which will decrease the perturbation if a > 1/3, while 
the third term will increase the perturbation. However, since a < 1, the exponential factor in this term will make 
this term insigniflcant in the future (a;3>l, ora^l). IfO<Q;<l/3, both terms will increase the perturbation and 

make it unstable. Eq.([7|) can be analytically solved by changing the variable — Mfec'^l-^"")^ : 



where C = (^V^^ e' 



' de 



d^ 



2c2 



and the solution of is 



3q: — 1 
2 - 2a 



ufe = 



(8) 



(9) 
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where v = \3a — l|/(2 — 2a), i^^ are modified Bessel functions and Ak, Bk are integration constants, which should 
be determined by the initial condition. When x ^ 1, namely a ^ 1, then ^ 1 and if a > 1/3 

then for a given mode, the perturbation will be a constant in the far future. Otherwise, if a < 1/3, <i>fc(a; — > oo) ~ 

oo, and if a = 1/3, ^k{x — *■ oo) ^ a; — > oo. Therefore, only if a > 1/3, the perturbation is stable. We 
also plot the evolution of perturbations in FiglT]and Fig[2]for an illustration. It should be noticed that B^ should 
depending on k, so for a given initial condition of Bk is different for each mode. 



3a 



a 



1/2 



(10) 




FIG. 1: RDE without matter: Evolution of perturbations with a = 0.3(left), 0.6(right), kH^ ^ — 0.5, 1.0, 1.5 and the initial 
condition $fco = 0.01, = 0-01. 




FIG. 2: RDE without matter: The most left figure shows the evolution of perturbations <l>fe with a — 0.5, 0.6, 0.7 and the initial 
condition $fco = 0.01, — 0.01. The middle and the most right figures show the evolution of with a — 0.6 and different 
initial conditions: $m = 0.01,0.02,0.03, = 0.01 in the middle one and $fco = 0.01, = -0.01,0.01,0.02 in the most 
right one. 

RDE with matter: If the matter component also exist in the universe, the background equation ([5]) becomes: 

+ ^7„oe-^^ (11) 

where h = H/Hq, flmo = Pmo/(3ffo) ^^'^ have used the energy conservation law for the matter. The solution to 
the above equation is 

- -^fl„oe-'^ + foe-('-i> , (12) 
Z — a 



,2 PR + Pra 



2h^ 



where /o can be determined by the condition h{0) — 



(13) 
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The energy density and pressure of RDE are 

PR 



2 - a 



"rnoe + Joe \ 



and 



Pi? = -Pn - = 



So, the equation of state at present is 



1 /2 



^.o = -3 --11(1- 



2 - a 1 — ilmO 



Again, only in the case of a < 1 , the universe will be accelerated expanding. Define 



Sp„ 



Spf 



2 - a 3f2mo 



(14) 



(15) 



(16) 



(17) 



where = 5pr/5{pr + pm) = p'r/^p'r + p'm) the total squared sound speed while c^j^ = SpR/Spu = p'r/ p'r is 
the squared sound speed for RDE only. By using the background equation (fTTj) and its solution (fT2|) . we obtain the 
following perturbation equation for RDE with matter: 



where 
and 



a \ 1 + 7 



l + f3 



$1. 



1 + 7 



f3 = r!,„oe-'"/j-' = 



2 - a flmo 



A 



3q: \ 1 + 7 



^0 ,(i-i) 



fc2i7-2e-2(i-i>$, = 0, (18) 



(19) 



A = 



2- a 



(20) 



In the limit of x ^ oo corresponding to the far future, then 7,/3 ^ and A ^ 0(1), therefore, eq.([T8]) comes back to 
eq.Q, which means the matter perturbation is not so important for the perturbation of RDE in the far future. Since 
we are interesting in the future behavior of the perturbations, we can make the following approximation: 
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and only keep linear terms of 7, /3, then the perturbation equation becomes: 



(3a -1)2 (2a-'^ + a-3)f}„oe"("~^)'^ 2- 



4a2 2a2(2a-l)/o 
where we have changed variable <i>fc to Wfc by 



3a /o 



3a»,„oe-(^-i)" ^ 
(2-a)(2a-l)/o 



fc2i/-2e-2(^-i)^ 



(21) 



Wfc=0, (22) 



= Mfc exp 



1 



2a 



27 

1 + 7 



Mfe exp 



3- 



(23) 



It seems that eg. (1221) can not be solved analytically, but we can solve it in the long- wave and short-wave limit as 
follows: For long-wave {kHQ^ <C 1) modes, ea. (l22]) reduces to 



(2a3 + a-3)r!„oe"(""')" (3a - 1)^ 



2a2(l-2a)/o 



4a2 



Uk = 0, 



(24) 
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whose solution is 

Uk = CifcJ^(a) + C2fcF^(a), 
where v = |3q; — l|/(2 — a), Ji/, Y^, are Bessel functions, Cik, C2k are integration constants and 



2a 



(2a3 + a- 3)»„o 
2a2(i_2a)/o 



1/2 



(25) 



(26) 



Therefore, when a; oo and a > 1/3, the perturbation <i>fc(x — » oo) « $c(fc), where $c(fc) is a constant for a given 

mode. Otherwise, if a = 1/3, '^k{x — * oo) ~ e~('^~°)^ ^ c>o, and oo) ~ a; ^ cx) if a = 1/3. Therefore, only 

if a > 1/3, the perturbation is stable in the long-wave limit. 
For short-wave {kH^^ ^ 1) modes, eq. ip^ reduces to 



3q!/o 



(27) 



and here we have neglected the second term proportional to e ^)^, since it is small in the far future(a; — > oo) and 
its solution is: 



so, Uk{x — > oo) ^ X. Here Difc, I?2fc are integration constants and 
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a 



2- a 



1 - a V 3a/o 



1/2 



(28) 



(29) 



Again, the perturbation is stable if a > 1/3 in the short-wave limit. If a = 1/3, ^k{x — > oo) ~ a; — + oo and if a < 1/3, 
^k{x oo) is exponential increased. We also numerically solve the eq. fTS]) and plot the evolution of (f>fe in Fig[3]and 
Figgl 




FIG. 3: RDE with matter: Evolution of perturbations $k with a = 0.3(left), 0.6(right), kHg ^ — 0.5, 1.0, 1.5 and the initial 
condition <l>feo = 0.01 and <E>J,o = 0.01. 

In conclusion, we have study the scalar perturbation of the Ricci dark energy (RDE) with and without matter in 
the flat FRW universe. To make the universe accelerated expanding, it requires the parameter a < 1 in RDE, and 
only if a > 1/3, the scalar perturbation is stable. Therefore, from theoretical aspect, the reasonable value for a is 
1/3 < a < 1. According to the joint analysis in ref.[l], the best-fit results with la uncertainly is a = 0-359^Q g25 and 
our result indeed gives a lower bound of a theoretically, which is more stringent than that from observations. 
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OjX) 020 



kn = 0-03 



FIG. 4: _RZ)i? with matter: The most left figure shows the evolution of perturbations $fc with a = 0.5,0.6,0.7 and the initial 
condition ^ko = 0.01, $5,0 ~ 0.01. The middle and the most right figures indicate the evolution of "1?^ with a = 0.6 and different 
initial conditions; $m = 0.01,0.02,0.03, $1,0 = 0.01 in the middle one and $ko = 0.01, = -0.01,0.01,0.02 in the most 
right one. 
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